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(54) Initializing model-based interpretations of digital radiographs 



(57) Automated computer aided diagnosis (CAD) 
processing of digital radiographs through model-based 
interpretation, with the initialization providing a set of in- 
itial parameters used by the model. The initial parame- 
ters can be selected based on expected pathology in 



the digital radiograph, and are optimized by the model 
to match features shown in the radiograph. The model 
can bean iterative model or a non-iterative model. Anal- 
ysis is performed on the interpretation result, so as to 
diagnose pathology shown in the radiograph. 
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Description 

[0001] The present invention relates to computer aid- 
ed diagnosis (CAD) of digital radiographs, particularly 
those that use deformable models in model-based in- 
terpretations of digital radiographs based on machine 
vision, and more particularly to initialization of such 
model-based interpretations. 

[0002] One promising advance in the field of compu- 
ter aided diagnosis (CAD) is the application of machine 
vision to digital radiographs. Of the various types of ma- 
chine vision techniques now available, CAD often em- 
ploys model-based interpretation based on deformable 
models of objects found in the radiographs. "Deformable 
models" are models that maintain the essential charac- 
teristics of the objects that they represent, such as bone 
and tissue structure, but deform to fit a range of exam- 
ples common to multiple different radiographs of differ- 
ent patients. 

[0003] One problem encountered with the use of mod- 
el-based interpretation methods is their initialization. 
Specifically, many model-based interpretation methods 
involve iterative searches in regions local to a current 
estimate for the model. At each iteration, the current es- 
timate is deformed slightly so as to provide a next iter- 
ated estimate. In general, the model converges through 
multiple iterations until it reaches a best estimate of the 
structures within the radiograph (with "best" meaning 
that further iterations are unlikely to lead to significantly 
different results). 

[0004] However, model-based interpretations, partic- 
ularly those that involve iterative searches, are prone to 
significant convergence errors if they have poor initiali- 
zation, for example, initialized from a bad starting loca- 
tion. Figures 1Aand 1B illustrate this situation. Figure 
1A illustrates a digital radiograph from a lateral lumbar 
spine examination, and further illustrates a model- 
based interpretation in the form of an active shape mod- 
el (ASM) attributable to Cootes and Taylor. See, for ex- 
ample, Cootes and Taylor, "Statistical model of appear- 
ance for medical image analysis and computer vision", 
Proceedings SPIE Medical Imaging, pp. 236-248 (Feb. 
2001). ASM employs an iterative search to deform a 
model from an initial position to a converged position. 
As seen in Figure 1 A, a poor initialization position leads 
to convergence at an incorrect estimate of the position 
of the lumbar vertebrae. On the other hand, as seen in 
Figure 1B, a good initialization position leads to correct 
convergence at accurate locations of the lumbar verte- 
brae. 

[0005] Good initialization parameters can be obtained 
by manual input by skilled radiologists or other medical 
personnel. However, manual initialization is counterpro- 
ductive to the goals of a fully automated CAD procedure. 
[0006] It is an object of the invention to address the 
foregoing situation through an automated analysis to 
obtain initial values for model parameters for a model- 
based interpretation of a digital radiograph. The inven- 



tion is based on the recognition of the inventors herein 
that all radiographs obtained in accordance with anyone 
specific radiographic protocol (for example, a lateral 
lumbar spine examination) will include distinctive re- 
5 gions common from patient-to-patient and from exami- 
nation-to-examination. 

[0007] Thus, according to one aspect, the invention 
obtains initial values for model parameters for a model- 
based interpretation of a digital radiograph obtained 

10 from a patient in accordance with a radiographic proto- 
col. A region of interest in the radiograph is identified, 
based on landmarks common to multiple different radi- 
ographs obtained with the same radiographic protocol. 
The region of interest is thereafter analyzed so as to cal- 

15 culate candidates for the initial model parameters. If 
needed, the candidates can be thereafter refined so as 
to disambiguate candidates with respect to repetitive 
structures found in the region of interest. 
[0008] In preferred embodiments, the landmarks 

20 comprise distinctive regions of high contrast within each 
different radiographic image produced by the same ra- 
diographic protocol. For example, in the case of a lateral 
lumbar spine protocol, radiographic images are charac- 
terized by a bright pelvic area, a dark area correspond- 

25 ing to a non-patient region beyond the patient's back, 
and a dark lung area. The dark lung area is separated 
from the dark non-patient region by a bright spine com- 
prising the region of interest. These landmarks together 
are identified using image enhancement techniques 

30 such as equalization, window leveling, and thresholding 
so as to define the region of interest. 
[0009] In further preferred embodiments, candidates 
for model parameters are calculated based on visually 
significant features in the region of interest together with 

35 spatial orientation of such features within the region of 
interest. Since the model might often contain repetitive 
structures which are difficult to disambiguate within the 
region of interest, the candidates are calculated merely 
to find some of the repetitive structures without distin- 

40 guishing one from the other. For example, in a lateral 
lumbar spine examination, the region of interest might 
consist of five vertebrae above the iliac bone. Initial pa- 
rameters for a deformable model might therefore consist 
of parameters that define five nearly-identical rectangu- 

45 lar regions, and it is therefore difficult to distinguish one 
of the vertebrae (and corresponding rectangular region) 
from another. Accordingly, the candidates often do not 
distinguish between the repetitive structures and often 
might contain a shift of one up or down from the actual 

50 position in the radiograph. 

[0010] Disambiguation is preferably performed rela- 
tive to the landmarks used to determine the region of 
interest, as well as relative to boundaries of the region 
of interest itself. For example, once significant vertebrae 

55 are identified in the region of interest for a lateral lumbar 
spine examination, and candidates calculated based on 
the significant vertebrae, the candidates can be disam- 
biguated by distance measurements relative to the dark 
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region of the lung and the bright region of the iliac bone. 
[0011] The initial values of the parameters may cor- 
respond to expected pathology in the radiograph. For 
example, in a situation where normal pathology is ex- 
pected, a "normal" initial model can be selected. Like- 
wise, in a situation where an abnormal pathology is ex- 
pected, an "abnormal" initial model can be selected. Al- 
ternatively, processing according to two or more of mul- 
tiple different sets of initial parameters can be conduct- 
ed, with automated selection of one of them being based 
on convergence of the deformable model. 
[0012] The model-based interpretation may be an in- 
terpretation based on deformable models, and can be 
an iterative model such as ASM described by Cootes 
and Taylor, and an active appearance model (AAM), al- 
so described by Cootes and Taylor. Alternatively, the 
model-based interpretation, might be a non-iterative 
model such as a model implemented through neural net- 
works or wavelet analysis of digital radiographic con- 
tent. 

[0013] Further preferred embodiments of the inven- 
tion involve automated CAD processing of a digital ra- 
diograph through calculation of initial values for model 
parameters for a model-based interpretation of the dig- 
ital radiograph, followed by revision of the initial values 
according to the model-based interpretation so as to ob- 
tain a best estimate of parameters that accurately model 
features found within the radiograph. Measurements 
are then obtained from the interpretation results so as 
to provide computer assisted diagnosis of pathology 
found in the radiograph. For example, in the case of a 
lateral lumbar spinal examination, CAD processing can 
be performed so as to diagnose kyphosis and lordosis, 
together with a quantification of the relative severity of 
these conditions. 

[0014] This brief summary has been provided so that 
the nature of the invention may be understood quickly. 
A more complete understanding of the invention can be 
obtained by reference to the following detailed descrip- 
tion of the preferred embodiment thereof, which is de- 
scribed by way of example only, in connection with the 
attached drawings in which: 

Figures 1 A and 1 B are representative digital radio- 
graphs showing the effects of poor initialization ver- 
sus good initialization; 

Figure 2 is a block diagram showing a teleradiolog- 
ical computer aided diagnostic (CAD) system; 

Figure 3 is a flow diagram showing CAD analysis 
according to an embodiment of the invention; 

Figures 4A and 4B are views for explaining desig- 
nation of training points in digital radiographs; 

Figure 5 is a view explaining automated search of 
image data in a digital radiograph so as to obtain 



accurate modeling of shapes therein; 

Figure 6 is a detailed flow diagram showing a meth- 
od for obtaining initial values for model parameters 
5 according to an embodiment of the invention; 

Figures 7 A through 7F are representative digital ra- 
diographs of different patients under the same ra- 
diological protocol; 

10 

Figures 8A through 8F are views for explaining im- 
age processing by which a region of interest is iden- 
tified; 

15 Figures 9A through 9F are views for explaining im- 
age processing by which candidates are calculated 
for initial model parameters; 

Figures 1 0A through 1 0F are views of the same dig- 
20 ital radiographs shown in Figures 7A through 7F, but 
with identified vertebrae mapped onto the image, 
for explaining disambiguation according to an em- 
bodiment of the invention; and 

25 Figure 1 1 is a flow diagram for explaining a second 
embodiment of the invention. 

[0015] Figure 2 is a generalized block diagram view 
of a teleradiological CAD (computer aided diagnosis) 

30 system. As shown in Figure 2, a teleradiological CAD 
system includes multiple hospitals and radiology cent- 
ers 10, 20 and 30, an administrative site 40, and a tel- 
eradiological CAD site 50, all interconnected through a 
wide area network 45 or over the network. A typical hos- 

35 pital includes digital radiography equipment 11 for ob- 
taining original digital radiographs from a patient, as well 
as a film scanner 1 2 for converting film x-rays into digital 
radiographic form. The hospital further includes PACS 
(picture archiving and communication system) worksta- 

40 tions 14 and 15, all of which intercommunicate with im- 
age database 16 over a network connection 17 which 
may be a local area network, a wide area network, or 
an intranet. A router 18 provides communications with 
other components of the teleradiological CAD system. 

45 [001 6] Hospital 20 and the additional hospitals and ra- 
diology centers 30 include similar architecture, although 
it is to be understood that these architectures are illus- 
trative only of the general nature of radiology centers. 
[001 7] Administrative site 40 administers the teleradi- 

50 ological CAD aspects of the system such as by accept- 
ing and routing requests to appropriate CAD sites as 
well as routing diagnostic information back to the re- 
questing site, as appropriate. 

[0018] Teleradiological CAD site 50 includes a CAD 
55 server 51 which communicates with PACS workstations 
52 and 53 as well as image database 55 over a network 
connection 57. Router 58 interconnects the teleradio- 
logical CAD site 50 to other components of the system. 
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[0019] In a representative operational aspect, hospi- 
tal 10 will obtain a digital radiograph using digital radi- 
ography equipment 11 or film scanner 12, which is 
stored on image database 16. A radiologist or other 
medical personnel using one of the PACS workstations 
14 or 15 issues a request for CAD services, which is 
handled by administrative site 40. Administrative site 40 
routes the request to teleradiological CAD site 50 where 
CAD server 51 services the request. Technicians at 
CAD site 50 are preferably involved with the CAD anal- 
ysis using one of PACS workstations 52 and 53. The 
image data itself might be obtained form image data- 
base 16, or might be transferred to image database 55. 
Preferably, the image is stored in DICOM format. The 
results of CAD analysis is routed back to hospital 10 
through administrative site 40, where it reaches the orig- 
inal requestor at one of PACS workstations 14 and 15. 
[0020] In this embodiment, the CAD system of the 
present invention resides on CAD server 51 . Of course, 
other embodiments are possible, such as embodiments 
where CAD is performed locally in PACS workstations 
14 or 1 5, or on the digital radiography equipment 1 1 it- 
self. However, the centralized approach of the present 
embodiment has advantages over distributed ap- 
proaches, such as the advantage of providing uniform 
diagnosis as well as the advantage of simplified update 
in the event of an upgrade in automated diagnostic ca- 
pabilities. 

[0021] Figure 3 is a generalized flow diagram showing 
CAD processing in CAD server 51 . CAD processing ac- 
cording to the invention is a model-based interpretation 
of digital radiographs in which a plurality of model pa- 
rameters define a mathematical model of features in the 
digital radiograph. Presently, two models are preferred, 
both attributable to Cootes and Taylor: the aforemen- 
tioned active shape model (ASM) in which shapes are 
modeled by matching grey-level intensity characteris- 
tics around model points so as to achieve a best match 
to corresponding points in training data; and active ap- 
pearance model (AAM) in which shape and texture are 
used to constrain object appearance during the search. 
Both ASM and AAM employ iterative searches so as to 
achieve a best match, and other iterative model-based 
interpretations can be employed. In addition, non-itera- 
tive model-based interpretations can also be employed, 
such as neural networks and the like. 
[0022] As shown in Figure 3, the initial values of the 
model parameters are calculated in step S301 using im- 
age data 61 representing the digital radiograph being 
subjected to CAD processing, as well as the model 63. 
The model is based on training data which is empirically 
derived data obtained beforehand, usually through anal- 
ysis of a large number of digital radiographs from the 
same radiographic protocol as that represented by im- 
age data 61. The empirical analysis is typically under- 
taken by skilled physicians by the means of hand anno- 
tations of training images. The physicians draw land- 
mark points in each of the training images, then addi- 



tional points along boundaries are generated by inter- 
polation between the landmark points. Figures 4A and 
4B illustrate how a training shape is generated by this 
method, using the example of ASM. The image in Figure 

5 4A shows landmark points that a skilled physician draws 
by hand, whereas the image shown in Figure 4B shows 
the shape resulting from interpolating along boundaries. 
[0023] Although annotation of hundreds of training im- 
agesby hand using skilled physicians is a time-consum- 

10 ing task, it is preferred to automatic or semi-automatic 
methods of annotation since the accuracy and reliability 
of skilled physicians yields improved training data. 
[0024] Reverting to Figure 3, Step S301 generates in- 
itial values for the model's parameters as explained 

15 more fully below in connection with Figure 6. The gen- 
erated initial values are shown diagrammatically at ref- 
erence numeral 62. After initialization in step S301, 
model parameters are optimized in step S302. The pre- 
cise method of optimization depends on the particular 

20 model-based interpretation being employed. In the case 
of ASM, an iterative search is performed so as to inspect 
changes in image intensities along profiles normal to the 
model boundary through each model point. For a given 
model point, the grey-level intensity (or its derivative) is 

25 sampled along a profile of k pixels on either side of the 
point in the image. This is illustrated in Figure 5, which 
shows line segments extending outwardly from "0" 
marks, and along which image intensity characteristics 
are sampled. A multivariate Gaussian model of grey-lev- 

30 el intensity (or derivative) samples is created for each 
model point using profiles from the training set. During 
the ASM search, each model point is moved along the 
profile to achieve a "best" match in grey-level intensity 
characteristics, with "best" meaning that no significant 

35 improvements in accuracy are obtained through further 
iterations. Intuitively, if the model boundary corresponds 
to an edge, the aforementioned search process will lo- 
cate the most similar edge along the profile. After up- 
dating all model point positions, new model parameters 

40 are found to fit the model shape to the new shape. This 
process, of moving model points to best match the im- 
aging characteristics and then updating the model 
shape parameters, repeats until convergence is 
achieved since no significant change in point positions 

45 is obtained with further iterations. 

[0025] Asimilar approach is undertaken using model- 
based interpretation according to AAM, although AAM 
uses not only the shape of features within an image but 
also uses the description of texture across the object. 

50 [0026] After model parameters have been optimized 
in step S302, step S303 operates to extract diagnostic 
information from the optimized parameters. For exam- 
ple, in the case of a radiographic protocol of a lateral 
lumbar spine examination, automatic measurements 

55 are made of disk space, vertebral height, etc. Extraction 
of diagnostic information such as by automatic meas- 
urement proceeds based on the interpretive model used 
as well as the radiographic protocol employed. For ex- 
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ample, in radiographic protocols involving the forearm, 
foot, or hand, infant hip, infant foot, and leg x-rays, au- 
tomatic measurements are made of significant image 
characteristics common to those protocols. As a further 
example, in the case of three-dimensional radiographic 
images, three-dimensional measurements may be 
made such as measurements that might differentiate 
between scoliosis, kyphosis and lordosis. Finally, in step 
S304, a diagnostic report is output. 
[0027] All model-based interpretations of imagery 
need good initialization of model parameters so as to 
ensure convergence. ASM starts with an initial shape, 
based on which imaging characteristics are extracted 
and used to improve the shape. AAM starts with some 
initial appearance model parameters, which provide 
both initial shape and initial texture. 
[0028] Figure 6 is a flow diagram showing initialization 
of model-based interpretations according to an embod- 
iment of the invention (step S301 in Figure 3). Briefly, 
Figure 6 illustrates a technique for obtaining initial val- 
ues for model parameters in a model-based interpreta- 
tion of a digital radiograph obtained from a patient in ac- 
cordance with a radiographic protocol, in which the mod- 
el-based interpretation revises the initial values to ob- 
tain revised values for the parameters based on content 
of the digital radiograph so as to model features therein. 
As shown in Figure 6, a region of interest is identified in 
the radiograph, wherein the region of interest is identi- 
fied based on landmarks common to multiple different 
radiographs obtained from the same radiographic pro- 
tocol. The region of interest is thereafter analyzed so as 
to calculate candidates for the initial model parameters. 
If the candidates are unambiguous, then the candidates 
are used as the initial model parameters for the flow 
shown in Figure 3. On the other hand, if the candidates 
are ambiguous, then the candidates are disambiguated 
with respect to repetitive structures found in the region 
of interest. 

[0029] Turning more specifically to Figure 6, the in- 
ventors herein have recognized that all radiographs ob- 
tained in accordance with any one specific radiographic 
protocol (for example, a lateral lumbar spine examina- 
tion), will include distinctive regions common from pa- 
tient-to-patient and from examination-to-examination. 
These distinctive regions form landmarks common to 
multiple different radiographs, and these landmarks can 
be used to identify a region of interest in any one par- 
ticular radiograph. For example, Figures 7A through 7F 
illustrate different radiographs from different patients all 
for the same protocol, which, in this case, is a lateral 
lumbar spine examination. Although the patients are dif- 
ferent and the examinations are different, each of Fig- 
ures 7A through 7F include common landmarks in the 
form of a bright pelvic area, a dark area corresponding 
to a non-patient region beyond the patient's back, and 
a dark lung area. These three areas surround the lumbar 
spine, which is the region of interest for this particular 
radiographic protocol. Identification of the region of in- 



terest based on these common landmarks is described 
in more detail in connection with steps S601 through 
S606. First, recognizing that original images are rarely 
perfect, window leveling is performed in step S601 . This 
5 is illustrated in Figure 8B in connection with an original 
digital radiograph shown in Figure 8A (which also cor- 
responds to the radiograph shown in Figure 7A). 
Thresholding is also applied in step S601 (shown in Fig- 
ure 8C) with the intent being to locate the iliac bone. The 
spine is detected in the top part of the image, adjacent 
the dark lung area (step S602). The right boundary of 
the region of interest is then detected by connecting the 
right boundary of the spine with the top of the iliac bone 
(step S603), as illustrated in Figure 8D. The image is 
rotated in step S604 in order to make the boundary ver- 
tical, which provides simplified processing based on 
horizontal and vertical lines and structures. A rotated im- 
age is shown in Figure 8E. A rectangular region of in- 
terest is then obtained based on empirically derived in- 
formation concerning the radiographic protocol. In the 
case of a lateral lumbar spine examination, it has been 
determined that suitable regions of interest are approx- 
imately 1/4 of the image width, and based on the lower 
level of the lung bottom, a rectangular region of interest 
is obtained for the lumbar vertebrae. Additional window 
leveling is performed within the region of interest to ac- 
centuate features therein (step S606). The leveled re- 
gion of interest in depicted in Figure 8F. 
[0030] After the region of interest is identified, steps 
S608 through S613 operate to analyze the region of in- 
terest so as to calculate candidates for the initial model 
parameters. In the case of a lateral lumbar spine exam- 
ination, the model parameters define the shape of five 
vertebrae within the region of interest. Since these ver- 
tebrae are in relatively fixed position to each other, the 
center of one vertebra in the image provides sufficient 
information to obtain candidates for the model parame- 
ters. Extra information about the vertebrae can further 
improve initialization of model parameters. 
[0031] Thus, in step S608, the image within the region 
of interest is sharpened, and a simple 3x3 filter mask is 
applied to locate horizontal lines (step S609). The re- 
sults of these steps are depicted in Figures 9A and 9B, 
respectively. A brightest horizontal line is found, and its 
width calculated (step S610), this brightest horizontal 
line being assumed to correspond to a vertebral edge. 
Then, two other neighboring edges are detected (step 
S611). The edge detected depends on the location of 
the brightest vertebra relative to the lung and the iliac 
bone. In particular, the neighboring edges are searched 
downward to the first edge, or upward to the first edge, 
or one up and one down depending on this location. In 
any of these cases, the original edge plus the two neigh- 
boring edges define one vertebra and the edge of a ver- 
tebra adjacent thereto. Results of edge detection is 
shown in Figure 9D. 

[0032] A fourth edge is then detected (step S612) so 
as to outline two adjacent vertebrae, as shown in Figure 
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9E. These edges are rotated back to the original orien- 
tation (step S61 3), as shown in Figure 9F. The center of 
one of the outlined vertebrae could be chosen to initial- 
ize the model parameters. 

[0033] With many images, the resulting approximate 
location is sufficient to define candidates for the initial 
model parameters. For example, when applied to a ra- 
diographic protocol involving cartilage of the knee, the 
shape defined by the initial model parameters is suffi- 
ciently distinct that the candidates are usable for the in- 
itial model parameter. However, in other circumstances, 
particularly one in which there are multiple repetitive 
structures that are similar to each other, disambiguation 
may be necessary to disambiguate the candidates with 
respect to these repetitive structures. 
[0034] That situation presents itself in lateral lumbar 
spinal examinations, where the five vertebrae are repet- 
itive rectangles that are often difficult to distinguish. Fig- 
ures 10A through 10F illustrate this situation, in which 
steps S601 through S613 have been applied to the six 
digital radiographs originally illustrated in Figures 7A 
through 7F. As seen in Figures 10A through 10F, steps 
S601 through S613 have resulted in candidates that 
identify different vertebrae in each patient. As a conse- 
quence, disambiguation is needed for this protocol. 
[0035] If disambiguation is needed, it is performed in 
step S61 5. In the case of lateral lumbar spinal examina- 
tions, disambiguation is performed by identifying each 
vertebra using its relative location to the lung and iliac 
bone. 

[0036] Figure 11 illustrates a second embodiment of 
the invention in which multiple different sets of initial pa- 
rameters are generated, one each for respectively dif- 
ferent pathologies in the radiographic protocol in ques- 
tions. Likewise, multiple different models are provided, 
one for each of the different sets of pathologies. For ex- 
ample, in the case of a lateral lumbar spinal examina- 
tion, different pathologies can include scoliosis, kypho- 
sis and lordosis. The model-based interpretation uses 
associated models and initialized model parameters so 
as to obtain a corresponding multitude of converged 
model parameters (step S1102). Each set is analyzed 
to determine which has converged the best, and con- 
verged model parameters for that set are chosen (step 
S1103). For the chosen convergence, diagnostic infor- 
mation is extracted (step S1104) and a diagnostic report 
is output (step S1105). 

[0037] The invention has been described with respect 
to particular illustrative embodiments. It is to be under- 
stood that the invention is not limited to the above-de- 
scribed embodiments and that various changes and 
modifications may be made by those of ordinary skill in 
the art without departing from the scope of the invention. 



Claims 

1. A method for obtaining initial model parameters for 



a model-based interpretation of a digital radiograph 
obtained from a patient in accordance with a radio- 
graphic protocol, wherein the model-based inter- 
pretation changes the model parameters based on 
5 content of the digital radiograph so as to model fea- 

tures therein, said method comprising: 

identifying a region of interest in the radiograph, 
wherein the region of interest is identified 
10 based on landmarks common to multiple differ- 

ent radiographs obtained with the same radio- 
graphic protocol; and 

analyzing the region of interest so as to calcu- 
late one or more candidates for initial model pa- 
15 rameters. 

2. A method according to Claim 1, wherein the land- 
marks comprise distinctive regions of high contrast 
within multiple different radiographic images pro- 

20 duced by the same radiographic protocol. 

3. A method according to Claim 2, wherein the radio- 
graphic protocol is a lateral lumbar spine protocol, 
and wherein the landmarks include a bright pelvic 

25 area, a dark area corresponding to a non-patient re- 
gion beyond the patient's back, and a dark lung ar- 
ea, which is separated from the dark non-patient ar- 
ea by a bright spine comprising the region of inter- 
est. 

30 

4. A method according to Claim 1 , wherein said iden- 
tifying step comprises image enhancement tech- 
niques including equalization, window leveling, and 
thresholding so as to define the region of interest. 

35 

5. A method according to Claim 1 , wherein the candi- 
dates for initial model parameters are calculated 
based on visually significant features in the region 
of interest together with spatial orientation of such 

40 features within the region of interest. 

6. A method according to Claim 5, wherein the radio- 
graphic protocol is a lateral lumbar spine examina- 
tion, and the initial parameters for the deformable 

45 model define five nearly-identical rectangular re- 
gions corresponding to five vertebrae above the il- 
iac bone, and the candidates for initial model pa- 
rameters are calculated without regard to which of 
the five rectangular regions corresponds to one of 
50 the vertebrae. 

7. A method according to Claim 1 , further comprising 
the step of disambiguating the candidates for initial 
model parameters with respect to repetitive struc- 

55 tures found in the region of interest. 

8. A method according to Claim 7, wherein disam- 
biguation is performed relative to the landmarks 
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used to determine the region of interest. 

9. A method according to Claim 8, wherein disam- 
biguation is also performed relative to boundaries 

of the region of interest itself. 5 

10. A method according to Claim 9, wherein the radio- 
graphic protocol is a lateral lumbar spine examina- 
tion, and the initial model parameters are disam- 
biguated by distance measurements relative to a 10 
dark region of the lung and a bright region of the 
iliac bone. 

11. A method according to Claim 1 , further comprising 

the step of selecting an initial set of model param- 1$ 
eters from among multiple different initial sets cor- 
responding to multiple different pathologies. 

12. A method according to Claim 1 , further comprising 

the steps of: 20 

selecting multiple different models correspond- 
ing to multiple different pathologies; 
obtaining an initial set of model parameters for 
each different model according to Claim 1 ; 25 
changing each initial set of model parameters 
according to the model-based interpretation; 
and 

selecting one set of model parameters based 

on convergence of all models in the model- 30 

based interpretation. 

13. A method according to Claim 1 , wherein the model- 
based interpretation is based on an iterative model. 

35 

14. A method according to Claim 1 , wherein the model- 
based interpretation is based on a non-iterative 
model. 

15. Automated CAD processing of a digital radiograph 40 
through identification of a set of initial model param- 
eters for a model-based interpretation of the digital 
radiograph according to any of Claims 1 to 14, com- 
prising: 

45 

changing model parameters according to the 
model-based interpretation so as to obtain a 
best estimate of features found within the radi- 
ograph; and 

analysis of the interpretation results so as to 50 
provide computer assisted diagnosis of pathol- 
ogy found in the radiograph. 
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Figure 1A Figure IB 
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